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Abstract 


Variational  expressions  for  the  scattering  matrix 
of  a  wire  grating  obtained  in  a  previous  report  (EM-18)  are 
evaluated  for  the  case  of  circularly  cylindrical  wires.  The 
single  scattering  approximation  is  used  as  a  trial  field. 
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1,  Introduction 

In  a  previous  report  ,  stationary  expressione  for  the  scattering 
matrix  elements  of  a  wire  grating  were  derived.  These  expressions  do  not  involve 
any  asstunptions  ahout  the  size  or  shape  of  the  wires.   In  that  report,  these 
stationary  expressions  were  evalviated  for  the  case  of  thin  wires,  in  which  the 
shape  does  not  play  any  role.   In  the  present  report,  which  should  he  considered 
as  an  appendix  to  the  previous  one*,  these  expressions  will  he  evaluated  for  the 
specific  shape  of  the  circular  cylinder,  hut  with  no  restrictions  on  size.  The 
trial  current  distrihutions  will  he  taken  from  the  solution  of  the  prohlem  of 

diffraction  of  a  plane  wave  hy  a  single  circular  cylindrical  wire.   The  paper  of 

2 
Schaefer  and  Reiche  derives  the  intensities  of  the  various  orders  from  that 

solution  directly.   In  this  paper  the  results  are  improved  hy  using  the  current 
distribution  obtained  hy  neglecting  interaction  in  a  variational  expression. 
The  Green's  function  to  he  used  in  this  analysis  is  identical  with 
that  derived  previously,  hut  it  is  expressed  differently.  Thus,  instead  of  look- 
ing for  a  solution  of 

00 

(V^  +  k^)G(s;x,z;x',z')  =  XI  exp(insiragin9)6(x-x>-np)6(z-2')   (?) 

-00 

with  the  appropriate  periodicity,  we  recognize  that  the  solution  of  the  wave 
equation  with  the  specified  set  of  inhomogeneous  terms  is  the  sum  of  the  solutions 
of  the  wave  equation  with  only  one  inhomogeneous  term.   Thus  we  have 

00 

G(s;r;r')  <=  Y.    ^^   (s;r;r')  (yUa) 

-00 


(V^+k^)G'^  =  exp(ins'<asine^)6(x-x'-na)6(z-z<)  (7^^) 

This  equation  is  readily  solved  in  polar  coordinates, 

_n   i  Tt(l)  /,  I-*   — ^   -*i  I  \    inskasin  9i  /-,n\ 

G  =  ^-  H^  '  (kjr  -  na  -  r'l  ;  •  9         ^     ,  (7^) 

2.  TE  Wave  Problem 

We  shall  now  take  up  the  TE  wave  problem.   First  let  us  calculate  the 
current  distribution  excited  in  a  single  cylinder  of  radius  R  by  a  pleuie  wave.   In 
polar  coordinates,  with  6=0  the  direction  of  the  incident  wave,  we  have 

The  equations  in  this  report  axe   niunbered  consecutively  with  those  in  HI-lS, 
References  to  equations  (1)  through  (73)  will  be  found  in  El^I-lg. 
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rt  =  0    +  0 

=  e^p  (ikr  cos  e)  +  f;  A  H^^\kx)e^''® 


(76) 


00 


where  the  scattered  field  has  "been  expanded  in  a  series  of  eigenfunctions,  the 


coefficients  A  to  "be  determined.  Expanding  the  incident  field  gives 
n 


00 

-00 


(1)/ 


f-Y.  [i\(i")  *\V^^n  « 


ine 


(77) 


Since  the  "bo-undary  condition  is  that  f  vanish  on  the  cylinder,   r  =  R,  we  ottain 


0  =  i^'j   (kR)  +  A  H^^^kR) 
n  n  n 


(78) 


Patting  the  values  of  the  coefficients  A  thus  obtained  into  (77),  differentiating 
with  respect  to  r,  and  finally  evaluating  this  derivative  at  r  =  R, 


M^.e)  =  E  -S:-  w  {<^W).  v^)}  e^^^ 


(79) 


-00  H^-*-^kR) 
n 


where  V  is  the  Wronskian  of  the  two  functions.  Since  the  Wronskian  of  these  two 
Bessel  functions  is  2/TTRi,  we  take  as  the  current  distribution 

,n-l 


00 


.,^\       •<r-   2i"  "     in© 

-00  ttRH^  HkR) 
n 


(SO) 


Consider  now  the  grating  problem.  Let  the  plane  wave  be  incident  at  an 
angle  ©, ;  the  proper  current  distribution  is  then  1(0-©^^).   In  the  conjugate  problem. 


Diffracted  Wave 


Incident  Wave 


Figure  1. 
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the  wave  is   incident   in  the  direction  rr  +  9^,   so  that   the  dletritrution  is 
l(©-6„-n).      The   incident  fields  in  the  tv/o  cases  are   exp[ikr  008(0-6-  )]    and 
exp[-ikr  cos    (9-6    )]    respectively. 

We  can  now  proceed  to  the  evaluation  of  integrals   in  Eq.    (19) 
The  integrals   in  the  numerator  are 

f^  =y^  exp[ikE  c08(©-e2-TT)]    l(G-e^)  RdO  =  ^^(ej^.e,  +  Tt)  (gla) 

fg  =  /^   exp[ikR  co3(9-ej_)]  l{e-Q^-v)  RdO  =  t^{e^   +  TT,  ©^)  (glh) 

Bat,  hy  a  simple  change  of  variables,  (gla)  hecomes 

il=J^  exp(ikR  cose)  I(e  +  e^  +  TT  -  e^)  RdG  =  f(92  -  6^)  (82) 

a  function  of  one  variable  only.      Substituting  (gO)  in   (g2),   and  interchanging 
the  order  of  integration  and  summation,   we  get 

f(o  .e.)  =  f:      ^}"'         ein(TT^2-ei)y^  exp(ikRcose)e^-«d6  = 

2   ^      rSjTTH^^^kR)  -^^ 


,,      00     .n  in(ep-9i+n) 


QD       J     (kR)  .(r,      o     ) 

f{e^-e^)  =  -Ui  J2  — e^^^V^i^ 


-00  H   (kR) 
n 

It  is  readily  seen  that  the  n-th  term  of  ^i^^o'^l^   ^®  identical  with  the  (-n)-th 
term  of  f2(e^-e  ),  so  that  the  two  integrals  are  identical. 

The  next  step  is  to  evaluate  the  Green's  function  integral.   Instead  of 
retaining  the  form  given  by  (75)  and  integrating  around  the  same  wire  twice,  we 
replace  r'  +  na  by  p  thus  shifting  the  integration  with  respect  to  P  to  the  nth 
cylinder.   The  argument  of  the  Hankel  function  in  (75)  is  simply  the  wave  number, 
k,  times  the  distance  between  a  point  on  the  central  cylinder  and  another  point  on 
the  n-th  cylinder. 

In  the  0-th  term  of  the  Green's  function  integral  we  still  perform  the 
integration  around  the  same  cylinder  twice.  Using  the  addition  theorem  for  Bessel 
functions  ,  we  have 


H^^^(kl?-p|  )  =  f;  H^^\kR)J  (kR)e^^®-®'^  •    (gk) 

o  *— ^   m      m 


m=-oo 


-  k  - 


where  0  and  0'  denote  the  auigular  coordinates  of  r  and  P  ,  respectively.  We 
have  then, 

•  j23H^^\icR)j^(kR)y^i(e-e2-'rT)e^®Ede/fi(e'-e^)e-^®'Rde« 

m 
=  j5:H^'^kE)J  (kR)e^"'^®2^-^lW'"l(e)e^Rde/fl(0)e-^"'^d© 


Ji;4"('=H)j„(u»e'-(«2™-^i)!i5; 


m-1 


Ui 


-m-1 


H  (kR) 
— m 


H  (kR) 


(85) 


identical  with  the  integral  f.  Thus  we  see  that  if  we  neglected  the  rest  of  the 
terms  in  the  Green's  function  integral,  the  scattered  wave  ajnplitude  would  he 
essentially  f.  We  should  note  that  this  is  nothing  "but  the  far  field  amplitude 
scattered  hy  a  single  cylinder  from  the  same  plane  wave  in  the  6p  direction.  Neg- 
lecting the  other  terms  is  equivsilent  to  assuming  no  interaction  "between  grating 
cylinders.   Since  our  trial  current  dlstrihution  was  taken  as  the  exact  solution 
of  the  "no  interaction"  prohlem,  the  variational  expression  with  a  single  term 
in  the  Green's  function  retained  is  the  exact  scattering  matrix  element  for  that 
prohlem.  The  rest  of  the  terms  will  yield  an  estimate  of  the  effect  of  interaction, 
'('e  shall  now  evaluate  these  correction  terms.  In  order  to  do  this  we  need  an 
appropriate  expansion  of  the  Hankel  function.  Using  the  addition  formula  for  Bessel 
functions  of  order  zero  (for  definition  of  angles  and  distances  see  Fig.  2) 


Figure  2.  Two  wires  of  the  grating. 
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problem.  The  rest  of  the  terms  will  yield  an  estimate  of  the  effect  of  inter- 
action. We  shall  now  evaluate  these  correction  terms.   In  order  to  do  this  we 
need  an  appropriate  expansion  of  the  Hankel  function.  Using  the  addition  formula 
for  Bessel  functions  of  order  zero  (for  definition  of  angles  and  distances  see 
Big.  2) 

0  |J,         p. 


"•■"   Ll  LX 


^   ^ 


Applying  the  addition  theorem? to  H  » 


1^ 


eW^^kp)  =  EH^^(kna)j^(kR)e^ 


Hence 


and  the  n-th  term  of  the  Green's  function  integral  becomes 
^-insicasineig^^^  =  \fj  I (6-62-11  )H^^\k|?-P|  )i(e'-e^)R2dede' 

'   i:  ^  J  («^)J.(^)V.^^^yr  l(9-e2-TT)e-^^«Rde/fl(e.-9,)ei'^'F.d6. 
^^  .^*^l^(^lH(->-l)  V^f^     .  iii^^L- .  H   (n'^)e-^^(^2r^^)*^'^l 


Combining  this  result  with  that  for  the  (-n)-th  cylinder, 

J  J, 
-t  ,  -1 
H    H. 


?n  -  Sn+  *  g 


^_=.ME  ^  .  ^  (kR)H^^„(kna)[(-lfe-^^^'^^^^^l+(-l)V^'^""^^^«l]  (85) 


Summing  with  respect  to  n, 
Ag  =  f:g  =-Uii:  ^^(kR)f:H^i^nka)[(-lfe-'^^^^''^"l*(-l)^e-^^^^'^«l]   ^89) 

Restricting  ourselves  in  the  rest  of  the  discussion  to  the  somewhat  simpler  case 
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)f  normal  incidence,  Q-^  ~  ^* 


^S  =  -^i  E  [(-1  A(-1)T   it     .  Li  (kR)    .  f  H  ^    (nka)        .  C89a) 


3,  Evaluation  of  the  Hankel  Function  Series, 

The  series  involving  the  ratio  of  the  Bessel  to  the  Hankel  function  is 
rapidly  convergent  for  reasonatle  values  of  parameters.  Hence  both  the  previously 
considered  series  for  f  suid  the  douhle  series  are  satisfactory  for  numerical  com- 
putation. The  series  ^H   (nka),  however,  is  not  only  slowly  convergent,  hat  also 

fails  to  converge  absolutely  or  uniformly  with  respect  to  ka.  For  large  .n,H^^.jj(nka) 
behaves  like  exp(inka)/v/n'.  Although  the  series  T^  erp(inka)/v'n  can  be  expressed  in 
terms  of  the  generalized  Zeta  function  of  Riemann  ,  the  latter,  unfortunately,  is 
not  tabulated.  The  procedure  used  here  will  be  to  use  two  series,  a  linear  com- 
bination of  which  can  be  made  to  have  tv/o  terms  in  the  asymptotic  expansion  identical 
with  those  of  the  Hankel  function.   Consider 

11=1 
and 


00 


r(n)     n 


yp(x)  =  Z        '  ^"^    x"  (91) 

2     ^1   r(n+l/2) 

Using  Stirling's  theorem,  we  readily  obtain  asymptotic  series  for  the  Gamma  function 
quotients 

r(n) 


.^n-^/2   [  1  ^     1.     11      .H   ,.,     ]  .  (92b) 

r(n--l/2)  ^"^       58?^ 

1/2 
So  that  we  can  approximate  exp(inka)/  n  '        by  • 

In  fact  even  the  first  term  of  the  series  for  ^y^-^yg^i   n  =  1,   is 
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I  [  r(|)/  r(2)    +      r(i)/  r(|)]  =  1.0073.  only  0.770  off. 

Now  let  us  return  to  the  evaluation  of  the  series, 

5 
Introducing  the  Beta  function,  and  then  using  its  integral  representation 

^1  n 

and  ^ 

(    ^   „-l/2  v^  '^^2^'^^^^  n  _  -1/2  ^^,1     .  n 

-1/2  -r-  /I     n-l/T      v-1/2  n,  -1/2/1  dz 


(l-xz)^Z? 


1+,'  -^ 


=  .-V^/j5|      log-g.  (9U) 

Both  of  these  functions  are  singular  at  x  =  1  (or  at  ka/2n  an  integer)  and  they 
"behave  near  the  singularity  like  (l-x)~  '  .  A  parametric  plot  of  -k   [y,  (e   )  + 
y  (e  )]  is  shown  on  Fig.  3,  for  the  range  0  <  0  <  rr. 

Although  more  complicated  procedares  will  give  more  accurate  results, 

00 

in  this  paper  y^  H^   (nka)  will  he  approximated  hy* 
n=l  / 

(SHkar^/V'!/"  *^'l^  "  ''''^'[yiCe"")  *  y^C."'')].  M.   (89.)  kecome..  then. 
AC~-2(^  )'/'(E  [(-lAC-Dle''  "  "  ■'""'  ii  •  ^(Hi)|  by^)n^^^'-'^)]  .   (9, 


Finally,  the  right  hand  member  of  Eq.  (ig)  "becomes 


imd: 


-Ui  Y:   [j^(kH)/H^(kR)]e'°"'2  (96 

-i  m  

2»       ^.[(-1)  *(-l)"Je-^/(+»)n/2-in02(J^/H^)(Jy/H^)(kR)[y^(e^^^)-Ky2(e^^^)] 


1+ilL- 


yeFlki:  [  J2  J^(kH)/H^(kR)  e^''®2] 


n 


•  The  approximating  function  and  the  Hankel  function  series  have  the  sajne  first 
term  in  the  asymptotic  expansion. 
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^nx 


Fig.    3.     A  parametric  plot   of  |  [y^(e^®)  +  72(6^®)] 
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At  this  point  we  may  note  that  the  quantity  of  particuler  interest  in  (96)  is 
the  correction  factor 

2:  [(-ir*(-l)^]e-^^^-^^^^/2-ine2  i,  .  ii(^)^^^(^lka)^y^(^i^aj-] 

°TE~  ^^ —J  -— ■ (97) 

n   "n 

since  the  transmission  coefficients  "based  on  no  interaction  theory  have  "been  treated 

"before,  Near  critical  points  k  a  =  2mi,  the  correction  factor  dominates  the  be- 

^  1/? 

havior  of  the  transmission  coefficients,  which  go  to  zero  like  (k-k  )  '  , 

c 

1+,   TM  Vave  Pro"blem 

Let  us  take  up  now  the  TM  case.   We  derive  the  trial  current  distribution. 
as  above,  except  that  the  application  of  the  boundary  condition  that  the  normal 
derivative  vanish  gives  us,  instead  of  (78), 

iV(kR)  +  A^[H^^\kE)]  '  =  0       .  (98) 

The  current  distri"bution  is  now  proportional  to  the  function  on  the  boundary 

1(9)  =  f;    ^^"'     e^-«  .  (99) 

n=-0o  TTkEH'(kR) 
n 

The  normalizing  integrals  turn  out  to  be 

00  J'(kR)    ,     s 

f  =  .Ui  Y:    -^ e^^^V^a^  (100) 

-00  H'(kR) 
n 

instead  of  (83). 

Proceeding  to  the  evaluation  of  the  Green's   f\inction  integral,  we  note 

that 

m 
so  that,  carrying  out  the  integration  as  in  (85),  we  obtain. 

m 
again  identical  with  f,  for  the  same  reason. 
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For  the  remaining  terms  've  have 
(G^),,,  =  J  ^'  IZ  i^^*\c^)j;  (kfi)H  ^^(kna)e-^^^^'^'*^'^«^^^^^«l     (103) 


80  that. 

In  the  case  of  normal  Incidence,  the  correction  factor  is  then 
c^;.l.t^^i^ „   ^   1 (105) 


/2^^  TE  ^  (kR)e^^^2] 

1—   n     n  — ' 


n   n 


5,  Discussion  of  Results 

A  sample  calculation  was  performed  for  kR  =  2,  ka  =  8;  this  is  the 
case  of  the  half  open  grating  so  often  used  in  infrared  spectroscopy,  with  only 
one  order  of  diffraction  at  an  angle  51.76  .  The  correction  factor  for  the  TE 
case  was  found  to  "be   I.O5,  thus  indicating  that  for  most  practical  cases  the 
no  interaction  assumption  is  perhaps  satisfactory.   It  is  unfortvmate  that  the 
numerical  results  "based  on  this  assumption  are  not  available.  Most  autl  ors  having 
arrived  at  the  series  (83), state  that  it  does  not  converge  rapidly,  and  make  further 
approximations.  As  a  result  the  experimenter  is  left  with  the  verv  crude  Zirchhoff 
approximation  which  predicts  zero  intensity  of  even  order  spectra  in  the  half  open 
grating.  For  the  case  of  cylindrical  grating  wires  with  kR  =  U,  ka  =  I6,  the 
intensity  of  the  second  order  wave  (also  at  an  angle  ^l.lG   )  was  estimated  and 
found  not  to  he  negligible  compared  to  first  order  spectrum  intensity  for  the 
case  kR  =  2,  ka  =  8.  Thus,  in  millimeter  wave  spectroscopy,  the  contamination 
due  to  second  order  spectra  ceumot  be  neglected. 
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